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ABSTRACT 

Context. We study the peculiar interacting galaxy system of VCC1249/M49 located in the core of the Virgo B subcluster. Owing 
to a recent interaction between the dwarf galaxy VCC1249 and the halo gas of the elliptical galaxy M49, neutral hydrogen has been 
displaced from the interstellar medium of this dwarf into the Virgo intracluster medium. Observations also reveal multiple compact 
star-forming regions (aka HII regions) that are embedded in this HI cloud, with a projected separation up to 13 kpc from VCC1249 
in the northwest direction. 

Aims. Motivated by recent near-ultraviolet (NUV) imaging from the GALEX Ultraviolet Virgo Cluster Survey (GUViCS) of the 
VCC1249/M49 system that shows significant ongoing/recent star formation in the compact regions, we aim to constrain the origin of 
these outlying HII regions with a multi-wavelength approach. 

Methods. Using deep optical (;/, g, i, z) imaging from the Next Generation Virgo Cluster Survey (NGVS) and new Ha imaging 
obtained at the San Pedro Martir observatory together with Keck long-slit spectroscopy, we characterize the star formation rates, ages, 
and metallicity of VCC1249 and its outlying compact regions. Moreover, we analyze the color and luminosity profile of the galaxy to 
investigate its recent interaction with M49. 

Results. Our new observations indicate that VCC 1 249 underwent a recent interaction with M49 in which both ram-pressure stripping 
and tidal interaction occured. The joint action of the two mechanisms led to the removal of the HI gas from the interstellar medium 
of VCC1249, while the gravitational tides triggered the stellar tail and counter-tail of VCC1249. Our stellar population synthesis 
analysis reveals that the star formation in this galaxy was truncated around 200 Myr ago and that the outlying HII regions were born 
in situ x 10 Myr ago out of pre-enriched gas removed from the dwarf galaxy. These observations also reveal that interactions between 
central and satellite galaxies similar to those between VCC1249/M49 may be an effective way of dispersing metals into the halos of 
massive galaxies. 
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1. Introduction 

The role of the environment in shaping the observed proper- 
ties of galaxies (morphology, star formation, color, gas content, 
etc.) during their evolution is the subject of an ongoing debate. 
Among the multiple processes that operate in dense environ- 
ments (see the review by Boselli & Gavazzi 2006), ram-pressure 
stripping (Gunn & Gott 1972) is often invoked as the princi- 
pal mechanism acting in clusters, especially on dwarf galaxies. 
Indeed, due to their shallower gravitational potential, the dy- 
namical pressure that originates from the fast motion of these 
galaxies through the hot and dense intra-cluster medium (ICM) 
can easily exceed the gravitational binding force in the galaxy 
and efficiently remove atomic gas from their interstellar medium 
(ISM). The resulting sudden suppression of the star formation 
in these stripped galaxies, then, leads to their temporary trans- 
formation into post star-burst (K+A; Poggianti et al. 2004), 
and subsequently into dE galaxies (Gavazzi et al. 2010). Ram- 
pressure stripping can therefore be the leading mechanism for 
the migration of dwarf galaxies from the blue cloud to the red 
sequence (Boselli et al. 2008 a,b). 

Stripped neutral hydrogen is frequently observed in the lo- 
cal Universe. For instance, extended HI tails are observed in the 
Virgo cluster (Chung et al. 2007, 2009) associated with late-type 
galaxies within 1 Mpc projected distance from M87. Because 
the stripped gas always points away from the cluster center, 
these tails are a clear signature of the ram-pressure mechanism. 
Another typical feature of a ram-pressure stripped tail is the ab- 
sence of a stellar counterpart. Moreover, galaxies that are in- 
falling in rich clusters occasionally show stripped gas with as- 
sociated star formation (Cortese et al. 2007), as evident from 
their ultraviolet (UV) (Smith et al. 201 1), Ha (Yagi et al. 2010), 
and X-ray emission (Sun et al. 2010). A remarkable example is 
the dwarf irregular (dlrr) galaxy VCC1217, which shows an ex- 
tended tail of bright knots and of diffuse emission in the UV light 
(Hester et al. 2010; Fumagalli et al. 201 1), which is associated to 
star formation triggered by the interaction with the intergalactic 
medium. This interpretation is also consistent with the findings 
of hydrodynamic simulations (Kapfereret al. 2009). 

In addition to hydrodynamic processes, tidal interactions 
(e.g. Spitzer & Baade 1951; Valluri & Jog 1990; Due & Renaud 
2011) can efficiently remove baryons from galaxies, in partic- 
ular from the outer and less bound regions. However, because 
of the high relative velocities among cluster galaxies, these en- 
counters occur on shorter time scales than in the field, resulting 
in a less severe disturbance of the galaxy stellar and gaseous 
disks. Moreover, galaxies in clusters grow through merging and 
accretion processes, consistent with hierarchical models and as 
indicated by the extensive web of low-surface brightness fila- 
ments and tidal features found by Janowiecki et al. (2010) in 
the intracluster light around luminous ellipticals in Virgo (M49, 
M87, M86, M89). Perhaps the most spectacular example of an 
interacting system in the Virgo cluster is that between NGC4438 
and M86, whose morphology and gas content require that tidal 
interaction and ram-pressure stripping are simultaneously occur- 
ing (Kenney et al. 2008). Similarly, the perturbations observed in 
NGC4654 and NGC4254 can be modeled with a close encounter 
with a nearby companion combined with ram-pressure stripping 
(Vollmer 2003; Vollmer et al. 2005). 

Environmental processes may also play a role in the forma- 
tion of ultra compact dwarfs (UCD; Drinkwater et al. 2000), 
a new class of objects in clusters whose stellar masses exceed 
those observed in the brightest globular clusters only by small 
amounts (Hasegan et al. 2005; Mieske et al. 2008). However, 



their origin is controversial. One hypothesis is that they derive 
from stripping of dE galaxies in the cluster potential or the po- 
tential of its brightest members ("threshing scenario" Bekki et 
al. 2003). Alternatively, a merger of young and massive star clus- 
ters that formed during an interaction between gas-rich galaxies 
(Kroupa 1998; Fellhauer & Kroupa 2002) could account for the 
origin of UCDs. 

Similar arguments apply to the rare class of compact ellipti- 
cal (cE) galaxies, such as M32, whose compact nature has been 
ascribed to tidal stripping (e.g., Faber 1973; Bekki et al. 2001; 
Huxoretal. 2011). 

In this paper, we focus on the interaction between the din- 
galaxy VCC1249 ( aka UGC 7636) and the massive ellipti- 
cal M49. The system VCC1249/M49 has been studied previ- 
ously, following the detection of an HI cloud displaced from the 
galaxy in the direction of M49 (Sancisi et al. 1987; Patterson & 
Thuan 1992; Henning et al. 1993). Subsequently, McNamara et 
al. (1994) found a trail of debris that does not coincide with the 
HI gas and is projected northward from the dwarf galaxy^. Lee 
et al. (1997, 2000) furthermore showed the presence of young 
HII regions embedded in the HI cloud. The conclusions of these 
studies is that both ram-pressure stripping and tidal interaction 
are necessary to explain the morphology and the deficiency of 
HIinVCC1249. 

Our attention was again caught by the VCC1249/M49 sys- 
tem owing to a recent deep NUV GALEX image, obtained as 
part of the GALEX Ultraviolet Virgo Cluster Survey (GUViCS) 
(Boselli et al. 2011), showing an extended UV feature stretch- 
ing from VCC1249 toward M49, with some relative maxima 
in correspondence with the detached HI emission. By combin- 
ing these new data to high-quality optical imaging data from the 
Next Generation Virgo Cluster Survey (NGVS) (Ferrarese et al. 
201 1) at CFHT and with Keck optical spectroscopy of the exter- 
nal HII regions embedded in the HI cloud, we are able to conduct 
the first multi wavelength analysis of the properties of VCC1249 
and of the outlying HII regions. In particular, following a proce- 
dure similar to Fumagalli et al. (201 1) for VCC1217, we adopt a 
spectral energy distribution (SED) fitting technique to constrain 
the age of external HII regions and the epoch at which the star 
formation in VCC1249 was suddenly truncated, presumably by 
the interaction with M49. We stress that the ram-pressure mech- 
anism we are referring to throughout this work is not caused by 
the ICM within the Virgo Cluster, but by the hot and dense halo 
of M49, as described in Mayer et al. (2006). 

The organization of the present paper is as follows: in § 2, 
we discuss the observations and data reduction for each band. 
The properties of VCC1249 and of the external star-forming re- 
gions are studied in § 3 and § 4, respectively. In § 5 we show the 
results of the SED fitting analysis for both VCC1249 and the HII 
regions and in § 6 we summarize our conclusions and present a 
panoramic view of the Virgo's B subcluster updated accordingly. 

Throughout this paper we assume a standard cosmology and 
a distance modulus of 31.14±0.05 mag for M49 corresponding 
to the distance of 16.9+0.3 Mpc (Mei at al. 2007), fully consis- 
tent with 16.7+0.5 Mpc from Blakeslee et al. (2009) and with 
17.0 Mpc from Gavazzi et al. (1999). Magnitudes are given in 
the AB system throughout the paper. 
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Fig. 1. RGB image of VCC1249 (bottom-left) and of M49 (center) obtained by combining the NGVS images in the u,g,z filters. 
The difference in color between the red giant elliptical M49 and the blue dlrr VCC1249 is apparent. The faint tidal tail of VCC1249 
is also visible in the northwest direction. 



Table 1. Log of the observations. 



2. Observations and data reduction 



Imaging 



Telescope 


Filter 


Seeing 


Date 


Exp. Time 






(arc sec) 


(yyyy-mm-dd) 


(sec) 


CFHT 


m-NGVS 


0.96 


2009-02-27 


6402 


CFHT 


g-NGVS 


0.71 


2009-02-22 


3170 


CFHT 


i-NGVS 


0.61 


2009-02-18 


2050 


CFHT 


z-NGVS 


0.82 


2009-02-22 


4400 


GALEX 


NUV 


5.3 


2010-03-19 


1631 


GALEX 


FUV 


4.2 


2004-04-17 


104 


SPM 


Ha 6603 


1.5 


2011-04-25 


3x600 


SPM 


r' 


1.5 


2011-04-25 


300 


Spectroscopy 


Telescope 


Instrument 


Slit/Fiber 


Date 


Exp. Time 






(arc sec) 


(yyyy-mm-dd) 


(sec) 


Keck 


LRIS 


1.0 


2011-01-3 


2x1200 


Keck 


LRIS 


1.0 


2011-01-4 


900 


Lick 


KAST 


2.0 


2011-04-29 


2x1800 


SDSS 


SDSS 


3.0 


DR7 


2040 



Since 2009 we have collected new deep observations of the sys- 
tem VCC1249/M49, covering a wide stretch of the electromag- 
netic spectrum, including UV (NUV), optical imaging (u,g,i,z 
+ Ha), and optical spectroscopy. A summary of these observa- 
tions is presented in Table Q] In this section, we present a brief 
description of these observations and of the data reduction. 



1 This feature and the similarity in shape of the HI cloud and galaxy 
(see Figure [6} are the two strongest signatures for the ram-pressure 
mechanism highlighted by McNamara et al. (1994). 
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Fig. 2. (a) RGB image of VCC1249 obtained combining the NGVS images in the u,g,z filters. The outlying HII regions studied 
in this work are highlighted. Contamination from the light associated with the halo of M49 is clearly visible, (b) RGB image of 
VCC1249 obtained after subtracting a model for the light distribution in M49. (c) NUV contours superposed on the RGB image of 
VCC1249, after subtracting M49. (d) Enlargement of the region enclosed in the dashed yellow box in panel (b). Faint blue structures 
are visible near the highlighted HII regions. 



Images and spectra were analyzed using the STSDAS and 
GALPHO10 routines in the IRAItj and FUNTOOL43 packages. 

2 Developed for IRAF - STSDAS mainly by W. Freudling, J. Salzer, 
and M.P. Haynes (see Haynes et al. 1999) and adapted by L. Cortese 
and S. Zibetti to handle Ha data. 



3 IRAF is the Image Analysis and Reduction Facility made avail- 
able to the astronomical community by the National Optical Astronomy 
Observatories, which are operated by AURA, Inc., under contract with 
the U.S. National Science Foundation. STSDAS is distributed by the 
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Table 2. Aperture photometry corrected for extinction in the Milky Way of the two galaxies and of the star-forming regions. The 
photometry of M49 is computed within an aperture of about 5.7 arcmin radius (or ~ 1.6R e ; Ferrarese et al. 2006). 



ID 


Distance 


u 


8 


i 


z 


NUV 


FUV 




(arcsec) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0) 


(8) 


VCC1249 





15.60 ± 0.03 


14.73 ± 0.03 


14.19 ± 0.03 


14.11 ± 0.03 


17.46 ±0.15 


18.68 ± 0.18 


M49 


334.0 


10.28 ± 0.03 


8.83 ± 0.03 


7.81 ±0.03 


7.49 ± 0.03 


13.91 ±0.15 


14.82 ±0.17 


CI 


151.4 


21.33 ± 0.04 


21.31 ± 0.04 


21.53 ± 0.04 


22.12 ±0.04 


20.77 ±0.16 


20.80 ± 0.29 


C2 


126.1 


21.31 ±0.04 


21.41 ±0.04 


21.83 ±0.04 


22.51 ±0.05 


21.52 ±0.16 


21.09 ±0.32 


C3 


119.4 


23.53 ± 0.06 


23.27 ± 0.06 


23.15 ±0.05 


23.00 ± 0.06 






C4 


137.9 


23.09 ± 0.05 


22.56 ± 0.06 


22.33 ± 0.04 


22.28 ± 0.04 


(23.04 ±0.15) 


(22.48 ± 0.40) 


C5 


156.4 


22.91 ± 0.05 


22.69 ± 0.05 


23.00 ± 0.05 


23.28 ± 0.05 






C6 


70.7 


21.92 ±0.04 


21.94 ±0.04 


22.31 ±0.04 


22.10 ±0.04 


21.37 ± 0.16 


21.17 ± 0.33 


C7 


55.1 


23.17 ±0.06 


22.87 ± 0.06 


21.75 ±0.04 


21.22 ±0.03 






C8 


61.1 


22.31 ±0.06 


22.15 ±0.06 


22.14 ±0.05 


21.69 ± 0.04 


21.43 ±0.16 


20.98 ±0.31 


C12 


5.8 


22.17 ±0.06 


21.48 ±0.07 


21.11 ±0.06 


21.10 ±0.05 






C17 


19.2 


23.94 ± 0.09 


22.81 ±0.16 


22.03 ± 0.06 


22.08 ± 0.06 







(1) Naming of the star-forming regions as in Lee et al. (1997), (2) projected distance from the center of VCC1249. The photometric uncertainties 

are the quadratic sum of the ZP error and the Poisson error. 



The whole photometry was corrected for Galactic extinc- 
tion assuming E(B - V) = 0.022 from Schlegel, Finkbeiner 
& Davis (1998) and the Milky Way extinction law such that 
Afuv = %3%E(B-V),A NUV = %.1AE(B-V),A U = 4.90£(fi-V), 
A g = 3.70E(B-V),Ai = 2.00E(B-V) and A, = IA5E(B- V). 
Similarly a 6 % correction was applied to Ha for extinction in 
the Milky Way. 

2.1. Optical imaging: NGVS 

The system VCC1249/M49 was observed in the optical bands 
(u,g,i,z) as part of the NGVS0 under excellent seeing condi- 
tions (< 1 arcsec, see Table [TJ. The exposure times for each 
filter are listed in Table [TJ Images were reduced and cali- 
brated using a dedicated data reduction pipeline (Ferrarese et 
al. 2011) designed explicitly to recover faint, diffuse surface 
brightness features, reaching a surface brightness limit of /i g = 
29 mag arcsec -2 (2cr above the mean sky). The mean global sky 
background around VCC1249 was estimated and subtracted us- 
ing the GALPHOT tasks MARKSKY and SKYFIT. The former 
allows us to mark several rectangular regions on which the sky 
background is estimated; the latter computes the mean sky back- 
ground in each region using a sigma clipping algorithm, deter- 
mines the average between the regions and subtracts this con- 
stant value from the image. The sky regions are selected to avoid 
the halo of M49 and other bright objects. 

Because VCC1249 lies at 12 /, 30" ! 01 i .O + 07 rf 55"'46 s .O 
(J2000), only 5.6 projected arcmin from M49 (see Fig{T), it is 
difficult to derive accurate photometry of this galaxy and its star- 
forming regions from the optical images, (in particular in the 
reddest filters) without first removing the contamination from 
the halo of the nearby giant elliptical. We therefore subtracted 
the best-fit model of the M49 light distribution from the sky- 
subtracted images, which we obtained for each filter (u,g,i,z) 



Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under NASA 
contract NAS 5-26555. 

4 Developed for DS9 by the High Energy Astrophysics Division of 
SAO. 

5 https://www.astrosci.ca/NGVS/ 
TheJSfext_Generation_Virgo_Cluster_Survey/Home.html 
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Fig. 3. RGB cutouts of the regions CI (left), C2 (center), C6 
(right) obtained combining the u,g,z NGVS images. At the 
GALEX resolution of > 4 arcsec, CI and C2 are unresolved in 
the UV. The size of each image is about 10" x 10". 



using in sequence the tasks ELLIPSE and BMODEL after mask- 
ing VCC1249, the star-forming regions, and the bright stars. 
Panels (a) and (b) of Figure [2] show the RGB image before and 
after the subtraction of M49. To assess the quality of this sub- 
traction, we analyzed in each subtracted image the background 
level in circular rings centered at the position of M49, finding no 
residual flux within the uncertainty of the sky level. 

After removing the contamination from M49, we performed 
aperture photometry on the individual star-forming regions that 
are marked in Figure |2] The flux was measured in circular re- 
gions of about 3-4 arcsec in radius, subtracting the local back- 
ground determined within concentric annular regions of 5 to 10 
arcsec in radius. Throughout this work, we assume the nomen- 
clature introduced by Lee et al (1997). Notice the presence of 
other fainter blue features near the listed HII regions (see panel 
(d) Figure |5J, likely low surface brightness star-forming fila- 
ments or knots related to the interaction studied. Owing to their 
faintness and uncertain flux, these will not be further analyzed. 
For VCC1249 we also measured the surface brightness profile 
in each band (after masking the foreground stars and the back- 
ground objects), using a modified version of the ELLIPSE task 
in IRAF that fits elliptical isophotes. For consistency we per- 
formed the ellipse-fit on the g-band image and used the same 
elliptical isophotes for the other bands. During the analysis of 
the g-band image we maintained as free parameters the ellipse 
center, the ellipticity, and the position angle, and incremented 
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the ellipse semi-major axis by 5 arcsec at each step of the fitting 
procedure. 

2.2. UV imaging: GUViCS 

The system of VCC1249 and M49 was observed by GALEX in 
the NUV (1750-2750 A) as part of the GUViCSQ survey with an 
exposure time of 1630 sec. The NUV image, shown in Figure |4] 
reveals patchy extended emission between VCC1249 and M49. 
Compact emission is also seen in correspondence to some of the 
external HII regions. A shorter 104 sec exposure is available in 
the far-ultraviolet (FUV) band (1350-1750 A) from the GALEX 
all-sky survey. Owing to the shallower exposure, only VCC1249 
and the brightest outlying HII regions (CI, C2, C6 and C8) are 
marginally detected in the smoothed data. Because of the lower 
resolution of GALEX (> 4 arcsec), regions C3, C5, and C7 are 
not detected in the UV images. Furthermore, some regions that 
are unresolved at the GALEX resolution (e.g. C2) appear to be 
blends of two or more distinct HII regions resolved in the optical 
imaging (FigureO. CI is resolved into three components, one of 
which might be a background galaxy based on its red colors. 



12:30:00.0 



12:29:40.0 



■ ■if'. 




WW 

-■ V 



Fig. 4. NUV image of VCC1249 and M49 obtained from 
GUViCS. Multiple HII regions are trailing in the northeast di- 
rection from VCC1249. 



Similarly to what was done for the optical imaging, we per- 
formed aperture photometry for the galaxy and the outlying HII 
regions. Owing to the lack of deep IR data, the internal dust ex- 
tinction for VCC1249 is unknown. The IRAS 60 and 100 mi- 
cron upper-limits (GoldMind3 imply that A(FUV) is lower than 
~ 2 mag (Cortese et al. 2008). Moreover, the observed UV- 
optical colors of this galaxy are typical of blue-sequence dwarf 
galaxies (Boselli et al. 2008a; Cortese & Hughes 2009), sug- 
gesting that the amount of dust absorption is likely modest (< 1 



6 http://galex.oamp.fr/guvics/ 

7 Gavazzi et al. (2003), Galaxy On Line Database Milano Network 



mag). We tested the robustness of our results against the adopted 
value of dust attenuation, finding that our analysis is not affected 
if A(FUV) < 2. Thus, for simplicity, we assume A(FUV)-Q 
throughout this paper. For the HII regions instead we estimated 
the dust obscuration using the Balmer decrement. Because the 
measured I(Ha)/I(Hf3) ratio is 2.66 and 2.93 for C2 and C6, 
respectively, consistent with the expected Balmer ratio of 2.86 
(case B for T - 10000 K; Osterbrock 1989), no correction was 
applied. 

2.3. Ha imaging 

A 3x600 sec narrow band Ha image was acquired with the 
SPM3 2.12m telescope in April 2011 (see Table [TJ-These new 
data were combined with a 4x300 sec exposure acquired with 
the same instrument in 2001 that is available from the Goldmine 
database. A detailed description of the data reduction procedures 
for Ha images can be found in Boselli & Gavazzi (2002), and it 
is only briefly summarized here. The ON-band frame was ob- 
tained using a narrow band (80 A) filter centered at 6603 A 
that overlaps with the Ha emission line. The OFF-band frame, 
needed to subtract the stellar continuum, was acquired in a 
shorter exposure using the Gunn r-band filter. After bias sub- 
traction and flat-fielding, the OFF band image was renormal- 
ized to match the ON-band exposure. Finally, the OFF-band 
was subtracted from the ON-band to obtain the image of the 
continuum-subtracted NET Ha emission. The resulting NET and 
ON-band frames are shown in Figure|5] Although the normaliza- 
tion factor was determined assuming that no Ha flux should be 
found associated to stars in the field, some residual flux from 
the brightest stars in the field appears in the NET frame due to 
slight differences in seeing between the OFF and ON observa- 
tions (see Figure|5j- The absolute flux calibration was performed 
with repeated exposures of the standard spectrophotometric stars 
Feige34 and Hz44. 



|12:30:00.0"| |l2:29:5Q.¥| 




Fig. 5. Ha image of VCC1249 and the outlying HII regions ob- 
tained at the SPM telescope. Left: Ha net emission line. Right: 
Ha plus stellar continuum frame HII regions C 1 and C2 are la- 
beled. Stars with apparent residual NET emission (see text) are 
labeled SI to S5. Note the lack of substantial emission from 
VCC1249. 



In the final image, two sources (CI and C2) with an Ha 
surface brightness higher than cr„„„ = 2.0 x 10~ 17 erg cm" 2 
s arcsec -2 were detected (5<x above the background). As pre- 
viously, we measured the flux in circular apertures obtaining 



( http://goldmine.mib.infn.it I 



8 Observatorio Astronomico Nacional en San Pedro Martir, Baja 
California, Mexico. 



Arrigoni Battaia et al.: Stripped gas as fuel for newly formed HII regions in the encounter between VCC1249 and M49 



7 



(8.68 ± 0.33) x 1(T 16 erg cm- 2 s" 1 



and (1.49 ± 0.04) x 10" 



erg cm s for CI and C2, respectively. Since the narrow band 
filter is broad enough to include the bracketing [Nil] lines, the 
listed Ha fluxes were corrected accordingly, using [Nil] /Ha = 
0.1 1 for C2, as measured on the spectrum, while for CI we used 
an average ratio of about 0.1, as derived from the spectra of C2 
and C6. 



2.4. Spectroscopy 

We observed the HII regions C20, C4, and C6 on UT 2011 
January, 3 and 4 with the LRIS spectrograph (Oke et al, 1995) 
at the Keck I telescope on Mauna Kea. During the first night we 
obtained 2 x 1200 sec exposures of C2 and C4, using a 1" slit 
and with the instrument configured with the 600/4000 grism, the 
D560 dichroic and the 600/7500 grating, tilted to ensure contin- 
uous spectral coverage between the blue and red arms. A sin- 
gle exposure of 900 sec was obtained during the second night, 
with a similar configuration but using the 400/3400 grism. The 
2D images were reduced with the LowRedux pipeline^! which 
calibrates, extracts and fluxes the data. CI was subsequently ob- 
served on the UT night 2011 April, 29 with the KAST spec- 
trograph at the Lick observatory. We obtained 2 x 1800 sec ex- 
posures with the 600/4310 grism, the 600/7500 grating and the 
D55 dichroic using a 2" slit. For regions C12 and C17 the spec- 
tra taken in 3" apertures were retrieved from SDSS. 



3. Properties of VCC1249 

VCC1249 is a dlrr (Nilson et al. 1973) low-surface brightness 
galaxy, located 5.6 arcmin to the southeast of the giant ellipti- 
cal galaxy M49 (Kumar & Thonnard 1983), the brightest mem- 
ber of Virgo Cluster B, approximately 4 degrees (1.19 Mpc pro- 
jected distance) south of M87 (Virgo Cluster A). VCC1249 has 
an irregular morphology: the central part of the galaxy consists 
of several bright knots connected by bridges (Lee et al. 1997) 
and from this structure a tidal tail departs toward M49 in the 
northwest direction and a counter tail in the southwest direction 
(Patterson & Thuan 1992 and McNamara et al. 1994; see Figure 
|2|. Based on its photometric properties (Table|2), VCC1249 is a 
dwarf galaxy that lies in the blue cloud of the color-magnitude 
diagram of the Virgo cluster galaxies. 

VCC1249 has a systemic velocity of 390+30 km s _I (SDSS) 
B while M49 has 1001 + 12 km s" 1 (Schechter 1980). An HI 
cloud withMtf/ = (6.9 + 0.4)xl0 7 M o and radial velocity of 469 
km s -1 (McNamara et al. 1994), in agreement with the previous 
measurement of 472±4 km s _1 (Sancisi et al. 1987), belonging 
to VCC1249 and displaced toward M49, is detected in 21 -cm 
(see Figure|6]l. Conversely, VCC1249 itself shows no significant 
emission in the hydrogen line to a limit of Mhi < 4.2 x 10 7 M 
(Oosterloo & Shostak 1984). Quantified in terms of the HI defi- 
ciency (defm), defined as the logarithmic difference between the 
upper limit HI mass and the expected value for a galaxy of the 
same morphological type and size (Haynes & Giovanelli 1984), 
VCC1249 has defw > 0.92 and is among the most deficient 



9 LR1 in Lee et al. 2000 

10 http://www.ucolick.org/xavier/LowRedux/index.html 

11 This is the redshift reported by SDSS for regions C12 and C17. 
Huchra (1992) found 276±78 km s" 1 . Using this value, the relative ve- 
locity between VCC1249 and M49 would increase and the ram pressure 
would appear more effective. In turn, this would lead to a more rapid 
removal of the gas (58+21 Myr), leaving our conclusions unchanged, 
however. 
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Fig. 7. Surface brightness profile of VCC1249, obtained using 
ELLIPSE, letting the center, the ellipticity, and the position an- 
gle as free parameters. In red, we show the inner and outer scale 
lengths computed for each band adopting a two-component ex- 
ponential profile (superimposed in blue), r = yja • b, where a 
and b are the semi-major and the semi-minor axis of the isopho- 
tal ellipses, respectively. 



galaxies in the local Universe. This is not particularly surpris- 
ing, since at the projected distance between M49 and VCC1249 
of about 30 kpc, the X-ray emission is dominated by the hot (10 7 
K Forman et al. 1985) gas belonging to the halo of M49 with a 
density of more than 10~ 3 cirT 3 (Fabian 1985). 

In Figure [7] we show the surface brightness profile of 
VCC1249 as a function of the radius r = Va • b, where a and 
b are the semi-major and the semi-minor axis of the isopho- 
tal ellipses, respectively. The profiles are well fitted by an ex- 
ponential law with two components (for 0" < r < 30" and 
30" < r < 80"). The first component is characteristic of a disk, 
while the second one accounts for the presence of the tidal tail 
and its counter-tail (Patterson & Thuan 1992; Lee et al. 1997), 
which appear as an excess of light above the inner exponential 
fit, resulting in a greater scale length in the outer part of the pro- 
file. 

The color profiles of VCC1249 are shown in Figure [8] The 
u - i and g - i indices are as blue as u — i = 1 .20 and g - i — 0.39 
near the center, become redder moving outward up to 20 arcsec, 
then they flatten, consistently with the results of Lee et al. (1997) 
and Patterson & Thuan (1996). However, the NUV - u color 
shows an opposite trend: from NUV - m=1.8 near the center to 
0.8 at 40 arcsec. Note, however, that the uncertainties in these 
colors are large. 

Inspecting our Ha imaging, we see a lack of significant Ha 
emission associated to VCC1249 and to its tidal tails. Using 
an upper limit of F^ cl249 < 3.80 x 10~ 14 erg cirT 2 s , we 
infer an upper limit to the star formation rate (SFR) of about 
< 0.01 M Q yr -1 , using the calibration of Kennicutt (1998). 
Because Ha traces massive star formation, the lack of signifi- 
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Fig. 6. NGVS g image of VCC1249 (bottom-left) and of M49 (top-center) on which the NUV contours (red) and the HI contours 
taken from McNamara et al. (1994) (yellow) are superposed (the coordinates are precessed from B1950 originally used, to J2000). 
Note that the peak of the HI cloud nearly coincides with the peak of the NUV emission at the position of the region C2 (LR1 in Lee 
etal. 2000). 

sive stars. Using the FUV emission instead, we infer an SFR of 
about 0.005 ± 0.001 M yr 1 following Kennicutt (1998). The 
presence of far-ultraviolet continuum indicates that star forma- 
tion ocurred > 100 Myr ago, which combined with the lack of 
Ha emission, suggests that the star formation in VCC1249 was 
truncated recently. 



4. Properties of the external star-forming regions 

4. 1 . Morphology and photometry 

A complex of star-forming regions extends from the galaxy in 
northwest direction, up to 2.6 arcmin (about 13 kpc) far from 
VCC1249, as clearly visible from both the GALEX images and 
the NGVS images (Figure Ekb). 

Figure [9] shows the NUV - u colors of these star-forming re- 
gions compared to VCC1249. It appears that all the structures 
outside VCC1249 are much bluer than the galaxy (NUV — u - 
1.86 and u — i = 1.41). C12 and C17, belonging to the cen- 
tral region of VCC1249, are unresolved in the UV images and 
appear redder in the u - i color than all other regions, having 
(u - 0ci2 = 1-06 and (u — Ocn = 1.91, respectively. Conversely, 
the other regions have a mean color of (NUV - u) - -0.29 and 
(u - i) = 0.21. Only region C7 is as red as the regions belong- 
ing to the galaxy in u - i, but visual inspection of the imag- 
ing suggests that this is probably due to contamination from a 
background red source (see Figure |2}- Furthermore, the lack of 
a clear dependence of these colors as a function of the projected 
distance from the galaxy suggests that these regions are nearly 
coeval, regardless of the galactocentric distance. 
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Fig. 8. Color profile of VCC1249, obtained from the surface 
brightness profiles in Figure [7] 



cant Ha emission sets an upper limit on the age of the last star- 
formation event at > 20 Myr, i.e. the typical life-time of mas- 
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Fig. 9.NUV-U color of the studied HII regions and of VCC 1 249 
itself as a function of the projected distance from VCC1249. All 
outlying HII regions are bluer than VCC 1249. 

Only regions CI and C2 display Ha emission at the depth 
of our imaging. Their luminosities are 2.78 x 10 37 erg s _1 and 
4.78 x 10 37 erg s _1 , respectively, consistent with the faint end of 
the HII luminosity function (Kennicutt et al. 1989). Adopting 
the Kennicutt (1998) conversion, the inferred SFR is 2.20 x 
1(T 4 M yr 1 (CI) and 3.77 x 1(T 4 M Q yr 1 (C2). Using the 
FUV emission instead, we compute an SFR of about (7 + 2) x 
1(T 4 M Q yr 1 and (6 + 2) x 1(T 4 M yr" 1 , for CI and C2, re- 
spectively. The SFR of the other regions can be estimated from 
the FUV luminosity to approximately 1CT 4 M yr -1 per region 
(see Table EE 

4.2. Radial velocities and line properties 

The Keck spectra of C2 and C6 are shown in Figure \TU\ where 
we highlight the detected emission lines. These two spectra have 
a mean signal-to-noise ratio (S/N) of 3.7 and 2.20 respectively. 
The lower S/N in the spectra for the fainter HII regions C4 and 
CI, 0.9 and 1.4, respectively, prevented a detailed study of their 
emission properties, but we can derive an estimate of the radial 
velocity from the marginally detected Ha lines. Since the mea- 
sured I(Ha)/I(Hf3) ratio is about 2.66 and 2.93 for C2 and C6, 
respectively, consistent with the expected Balmer ratio of 2.86 
(case B for T = 10000 K; Osterbrock 1989) within the uncer- 
tainty, no reddening correction was applied to the spectra. 

We computed the radial velocities of C2 and C6 by averaging 
the redshift measurements obtained from the individual spectral 
lines. C2 and C6 have a radial velocity of 561±34 km s and 
533±53 km s _1 , respectively. Our redshift measurement for C2 
agrees with the value 577 ± 91 km s _1 from Lee et al. (2000). 
The radial velocities listed in Table [3] systematically increase to- 
ward M49 (1001 + 12 km s ). HI interferometric observations 
carried out with the VLA have shown that regions C1-C5 are 
apparently associated with an HI cloud of 6.9xl0 7 Mo at a 
recessional velocity of 469 km s _1 extending over 20 km s _I 
(McNamara et al. 1994). The VLA channel maps show a veloc- 
ity gradient across the HI cloud: the material in the northwest 
part has higher velocities (10-15 km s _1 ) compared to the center 
(McNamara et al. 1994). Having a bandpass of about 330 km s _I 
centered at 470 km s _1 , these observations (as well as the WSRT 
observation by Henning et al. 1993) cannot be used to check if 

12 The Kennicutt (1998) calibration applies to continuous star forma- 
tion over 10 8 yr or longer. Therefore, the estimated SFR of the HII re- 
gions should be approached with some caution. 

13 Obviously, the S/N turns out to be much higher in correspondence 
of the lines, exceeding 10 for the C2 and C6 regions. 



the association on the plane of the sky is also present in the ve- 
locity space. Lower resolution WSRT observations (Sancisi et 
al. 1987) and Arecibo observations (Hoffman et al. 1987) cover- 
ing a broader range in the velocity space, however, do not detect 
any 21 cm emission with an rms of 0.7 mJy at the redshift of CI 
or C4, and probably C2 and C6 (see Table 3). 

Table 3. Radial velocities for CI, C2, C4, C6, C12, and C17. 



Region 


Proj.Dist. 


Velocity 




(kpc) 


(km s" 1 ) 


(VCC 1249) 





390±30 


C12 


0.5 


390±30 


C17 


1.5 


390+30 


C6 


5.6 


533±53 


(HI) 


10.0 


469±3 


C2 


10.0 


561±34 


C4 


11.0 


656±73 


CI 


12.0 


716±106 


(M49) 


26.5 


1001+12 



The observed equivalent widths (EW) of the lines measured 
in C2 and C6 (adopting the convention that negative EW means 
emission) are listed in Table|5] From the ratios [OIII] A5007/H/3 
and [Nil] /16584/Ha', we determine that, as expected, C2 and C6 
lie on the HII regions arm of the BPT diagram (Baldwin et al. 
1981). No clear features related to shocks are seen in the spec- 
tra: no strong [OI] line is detected and all diagnostic ratios (in- 
volving [OI], [OIII], H/3, Ha, [SII], [Nil] lines) agree with the 
photo-ionization models that place C2 and C6 among HII re- 
gions (see e.g. Dopita & Evans 1986). 

Conversely, the SDSS spectra of C12 and C17 belong- 
ing to VCC 1249 exhibit no emission lines, but clearly show 
Balmer lines in absorption. In Table |4] we report the line equiv- 
alent widths. Consistent with our analysis of the star formation 
rate properties, the observed values are typical of k+a galaxies 
(EW(H5)>3 A, Poggianti et al. 2004, Dressier et al. 1999), post 
star-burst (PSB) galaxies whose star formation has been sud- 
denly truncated. These characteristics are not only found in C12 
and C17, but are representative of the entire galaxy. Indeed, the 
spectrum of VCC1249 (Patterson & Thuan 1992; Huchra 1992) 
exhibits no emission lines and strong Hf3 in absorption. 

Table 4. The Balmer absorption lines of C12 and C17. 





Wavelength 


EW (CI 2) 


EW (C17) 




(A) 


(A) 


(A) 


HS 


4101.7 


5 


4 


Hy 


4339.8 


6 


2 


W 


4861.0 


5 


4 


Ha 


6562.8 


1 


7 



4.2.1. Metallicity 

Since the [OIII] -14363 line is undetected in the spectra of C2 
and C6, we determined the oxygen abundances using the strong 
line methods. Several such methods, based on different line ra- 
tios with empirical, theoretical or "combined" calibrations, can 
be found in the literature. 

Since the abundances derived using these indirect methods 
are affected by up to ~ 0.7 dex uncertainties (Kewley & Ellison 
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Fig. 10. Spectrum of the HII regions C2 (top) and C6 (bottom) obtained at Keck. Dashed lines highlight the position of the lines in 
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2008), we calculated the chemical composition of the HII re- 
gions C2 and C6 using multiple ratios of strong lines3 For the 
R23 method, the observed log([NII]/Ha) > -1.1 puts C2 and 
C6 in the upper branch. 



14 fl 23 = {[Oil] 13727 + [0///] 114959, 5007) //^A4861, O n = 
{[OIII] 14959 + [OIII] 15007) / [Oil] 13727, [Nil] 16584/ ff a 16563, 
[OIII] 15007/ [2V7Z] 16584, [N//] 16584/ [0/7] 13727, and 
[#//] 16584/ [5//] 16720 



Then, following the procedure of Kewley & Ellison (2008), 
we homogenized the derived values to the Pettini & Pagel (2004) 
calibration to remove any systematic discrepancy between the 
various calibrations. The abundances derived using several cal- 
ibrations, before and after this conversion, are listed in Table [6] 
The inferred oxygen abundances show a peak to peak variation 
of ~ 0.2 - 0.3 dex. By averaging these different values, we ob- 
tained our best estimate of the metallicity 12 + log^) = 8.38 

for C2 and 12 + log(g) = 8.29 forC6. 



Arrigoni Battaia et al.: Stripped gas as fuel for newly formed HII regions in the encounter between VCC1249 and M49 
Table 5. Summary of the emission lines detected in the HII regions C2 and C6 



11 









C2 






C6 






A (rest) 


A (observed) 


Log Continuum 


EW 


A (observed) 


Log Continuum 


EW 




(A) 


(A) 


(erg cirr 2 s~ ] A -1 ) 


(A) 


(A) 


(erg cirr 2 s~ ] A -1 ) 


(A) 


[0//] 


3727.55 


3735.2 


-16.93 


-148.9 ± 1.1 


3734.8 


-17.04 


-29.5 ± 0.8 


HI 


3798 


3805.9 


-16.95 


-1.7 ±0.4 


- 


- 


- 


HI+Hell 


3834 


3842.4 


-16.95 


-1.4 ±0.5 


- 


- 


- 


[Nelll] 


3868.7 


3876.2 


-16.95 


-3.9 ± 0.4 


- 


- 


- 


H£+Hel 


3888.9 


3896.4 


-16.96 


-7.3 ± 0.4 


- 


- 


- 


He + [Nelll] 


3970.1 


3977.6 


-16.97 


-6.3 ± 0.3 


- 


- 


- 


[Sim 


4070 


4082.9 


-16.93 


-0.7 ± 0.3 


- 


- 


- 


[SII]2 


4078 


4087.1 


-16.93 


-0.9 ± 0.3 


- 


- 


- 


HS 


4101.7 


4109.6 


-16.95 


-9.8 ± 0.5 


- 


- 


- 


Hr 


4339.8 


4348.7 


-16.98 


-19.8 ±0.6 


4348.3 


-17.21 


-3.4 ± 0.5 


Hel 


4471.5 


4480.0 


-17.02 


-1.5 ±0.4 


- 


- 


- 


H/J 


4861.0 


4870.2 


-17.06 


-57.2 ± 0.6 


4869.4 


-17.43 


-13.6 ±0.9 


[<9///]l 


4958.9 


4967.9 


-17.06 


-24.5 ± 0.6 


4967.3 


-17.49 


-10.6 ± 1.1 


[0III]2 


5006.8 


5015.8 


-17.07 


-62.4 ± 0.8 


5015.2 


-17.52 


-34.3 ± 1.2 


NI 


5199 


5207.4 


-17.17 


-2.4 ± 0.6 


- 


- 


- 


Hel 


5875.6 


5886.3 


-17.23 


-6.6 ± 0.7 


5886.5 


-17.82 


-5.9 ±2.7 


[01] 


6300.3 


6311.7 


-17.25 


-5.9 ± 0.6 


- 


- 


- 


[NII]1 


6548.1 


6559.8 


-17.28 


-9.7 ± 0.6 


- 


- 


- 


Ha 


6562.3 


6574.8 


-17.28 


-252.9 ± 1.0 


6574.1 


-17.89 


-115.2 ± 1.5 


[NII]2 


6583.4 


6595.5 


-17.28 


-28.8 ± 0.5 


6595.3 


-17.89 


-10.8 ± 2.1 


Hel 


6678 


6690.5 


-17.25 


-3.2 ± 0.5 








IS 11] 


6717 


6728.8 


-17.25 


-38.8 ± 0.6 


6727.9 


-18.00 


-11.4 ±2.7 


ism 


6731.3 


6743.1 


-17.25 


-25.7 ± 0.6 


6742.6 


-18.05 


-17.1 ± 1.9 


Hel 


7065 


7078.4 


-17.26 


-0.9 ± 0.4 








[Arlll] 


7135.8 


7149.1 


-17.28 


-6.3 ± 0.5 








[01 l]l 


7319.5 


7333.5 


-17.27 


-4.0 ± 0.4 








[0/712 


7330.5 


7344.1 


-17.27 


-3.0 ± 0.6 









Table 6. Oxygen abundances of the HII regions C2 and C6. The 
values converted to the Pettini & Pagel (2004) calibration are 
given in parenthesis. 



Ratio 


C2 

12 + log (g) 


C6 

12 + %(g) 


Calibration 




8.69 (8.49) 


8.41 (8.34) 


McGaugh (1991) 


INII] 

loni) 


8.45 (8.24) 




Kewley & Dopita (2002) 


R23, O32 


8.85 (8.48) 


8.50 (8.19) 


Kobulnicky & Kewley (2004) 


R23 


8.81 (8.49) 


8.74 (8.42) 


Zaritsky et al. (1994) 


{Nil] 

H„ 


8.31 


8.25 


Pettini & Pagel (2004) 


INII] 


8.43 (8.31) 


8.35 (8.25) 


Denicolo et al. (2002) 



Assuming a solar abundance 12 + log (§) = 8.69 (Asplund 
et al. 2009), these values imply a subsolar metallicity in these 
two regions (Z = O.49Z for C2 and Z = 0.40Z Q for C6). 
Remarkably, this value is consistent with the metallicity derived 
for VCC1249 itself using the mass metallicity relation (see sub- 
section 15.11 ). which strengthens the hypothesis that the external 
HII regions were born in situ from the pre-enriched gas that has 
been stripped from VCC1249. 



5. SED fitting 

To constrain the star formation history of VCC1249 and its 
HII regions, we used the SED fitting technique presented in 
Fumagalli et al. (201 1). First, using the spectrophotometric evo- 
lution code PEGASE2.CPI we created a grid of synthetic spec- 
tra evaluated at different times for multiple input star forma- 
tion histories SFR(t). Then, using the SED-fitting code GOSSIP 
(Franzetti et al. 2008), we compared the observed photomet- 
ric points (and spectra if available) with the synthetic spectra. 
During this procedure, we analyzed VCC1249, M49, and the HII 
regions separately. 



5.1. The galaxy (VCC1249) 

As in Fumagalli et al. (201 1), the evolution of the spectral energy 
distribution of the galaxy was modeled assuming a Salpeter IMF, 
with a lower mass end of 0.1 M and an upper mass of 120 M Q . 
Initially, we assumed zero metallicity, then the ejecta of massive 
stars (Woosley & Weaver 1995), implemented in the PEGASE 
code, enrich the Furthermore, we assumed a delayed ex- 

ponential star formation history, dubbed "a la Sandage" (see 
equation 3 in Gavazzi et al. 2002, Sandage 1986), truncated 
as in Fumagalli et al. (2011) to simulate the stripping event on 



15 Projet d'Etude des GAlaxies par Synthese Evolutive. 

16 PEGASE uses the "Padova" stellar tracks, improved with 
the AGB phase (Groenewegen & de Jong 1993) and with 
the helium white dwarfs (Althaus & Benvenuto 1997). See 
http://www2.iap.fr/users/fioc/PEGASE.html for more details. 
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VCC1249 due to the interaction with M49. Analytically, this is 



SFR(t) 



£e 7 if t < t, rmc 
if t > t trunc 



where t trunc is the time from the onset until the end of star for- 
mation activity. 

In Figure QT| (adapted from Fumagalli et al. 201 1), we show 
an example of star formation history from our model library and 
we highlight the various timescales that are relevant for this anal- 
ysis. 



t 

age 




trunc 






t -t, 

aee trunc 






Today 






i-.-.-i i i 



5 10 
Timc(Cyr) 



Fig. 11. Example of star formation history from our library, with 
t=4 Gyr, t age - 8 Gyr and t trunc = 5 Gyr (adapted from Fumagalli 
etal. 2011). 



In this grid of synthetic spectra, the characteristic time scale 
of the burst t ranges from 1 to 20 Gyr with steps of 1 Gyr . The 
age of truncation t tnmc varies from 1 to 13 Gyr with steps of 1 
Gyr, while t age spans from to 13.5 Gyr with steps of 100 Myr. 
Our final library includes 35K spectra. During the fit, we did not 
fix the age for the onset of star formation activity and, as previ- 
ously discussed, we did not include the (unknown) dust correc- 
tion, since extinction is expected to be negligible for VCC1249 
(see figure 8 of Cortese et al. 2008 and 3-4 of Lee et al. 2009). 

With this library, we ran the GOSSIP software to evaluate 
the parameters that best reproduce the observed photometry and 
their associated probability distribution functions (PDFs). For 
VCC1249, we found that t age , t trunc , and t are not well con- 
strained individually and we were only able to constrain the 
two latter parameters to be < 2 Gyr. Conversely, the parame- 



ter t age - ttnmc, which represents the lookback time at which 
the truncation of star formation occurred, is better constrained 
tage - trmnc = 200^Jy[jMyr. Figure Q~2]shows the associated PDF. 
This best-fit value agrees with the PSB spectral features ob- 
served in VCC1249 that are expected to arise if the star forma- 
tion activity halted abruptly in the past 0.5 - 1.5 Gyr (Couch 
& Sharpies 1987). Furthermore, this is consistent with the ob- 
served lack of Ha emission and with the presence of faint FUV 
emission. 

Having inferred the star formation history of this galaxy, we 
can estimate its stellar mass by integrating the best-fit star for- 
mation history M star = (1.20 ± 0.16) x 1O 9 M . This value is 
consistent with what was reported by Lee et al. 2003, M star = 



0.001 
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Fig. 12. Probability distribution functions of the t age - t mmc pa- 
rameter for VCC1249. The truncation of star formation appears 
to be well constrained at about 200 Myr ago. 



o 

ac 
u 
o 

2 0.006 

o 



-h 0.003 



1 1 1 1 1 1 1 1 


ufcLAi M49 - 






— 1 — — 1 — 1 — 1 — — 1 — 1 — 1 — 


— i — i — i — — i — i — i — i — i — i — i — — i — i — i — 







?M(A) 



Fig. 13. Model SED for M49 (top panel) together with the ob- 
served NIR photometric points and an optical spectrum (green) 
from GoldMine. The best-fit model is for a Sandage star forma- 
tion history with t age = 12 Gyr and t = 2.8 Gyr. For VCC1249 
instead (shown in the bottom panel), we find that a truncated 



Sandage star formation history with a t, 
provides a good description of the observed SED. 



age 



of ~ 200 Myr 



9.55 x 1O 8 M . For a 10 9 M galaxy, the mass metallicity relation 
(Kewley & Ellison 2008 using the Pettini & Pagel 2004 calibra- 
tion) predicts 12 + log(^) ~ 8.35, which agrees well with the 
metallicity derived for C2 and C6. Again, this confirms a sce- 
nario in which the HII regions were bom from the gas stripped 
from VCC1249 and pre-enriched by previous episodes of star 
formation. 

Figure[13]shows the best-fit model and the SED of VCC1249 
(bottom panel) and of M49 (top panel) ( Xycci249 ~ 0.29; 
^M49 = 0-90). A detailed analysis of the elliptical galaxy is be- 
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Table 7. Derived parameters of the HII regions. The SFR is computed from Ha luminosity (left) and from FUV (right) for CI and 
C2 and only from FUV for the remaining regions using Kennicutt (1998) calibrations. Mass and age are obtained from the SED 
fitting analysis. 
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yond the scope of this study, but provides a test for our fitting 
tools. In this case, we complemented the NGVS measurements 
with J,H,K,U,B,V photometry and an optical spectrum taken 
from GoldMine. The best fit performed adopting a star formation 
history "a la Sandage" and a Salpeter IMF provides t age = 12+j 
Gyr and t = 2.8 ± 1 Gyr. Our determination agrees with the value 
of t age = 11.9 Gyr found by Idiart et al. (2007) and marginally 
agrees with t age -9.6+ 1 .4 Gyr of Thomas et al. (2005). 

5.2. The star-forming regions 

For the analysis of the external HII regions we performed again a 
modeling similar to the one presented in Fumagalli et al. (201 1). 
First, we computed a set of synthetic spectra with PEGASE2.0, 
assuming a single stellar burst with Salpeter IMF and slightly 
subsolar initial metallicity. As justified by the observed Balmer 
decrement, we neglected internal dust absorption. For C12 and 
C17, which belong to VCC1249, we performed the SED fitting 
using both the single-burst model and an exponentially decreas- 
ing SFR with t ranging from 10 to 2000 Myr. 

Table Q contains the best-fit parameters for the observed 
SEDs. The outlying HII regions are well described (e.g. Xci ~ 
0.63; Xce ~ 0.24) by a single burst of star formation, with age 
t < 30 Myr. Instead, the star formation history of C12 and C17 
is better described (x^ l2 = ^-^'<Xcn = 0.87) by an exponential 
with r < 20 Myr and with an age about a factor of 100 higher 
than that of the HII regions outside the galaxy, such as C2 or 
Cl. This result is consistent with the mean ages of » 10 7 yr and 
10 8 yr obtained by Lee et al. (1997) assuming the (B - V)-age 
relation given by Bruzual & Chariot (1993). 

Figure [14] shows the SEDs of the regions, along with the 
PDFs of the age parameter that appears to be well constrained. 
As previously inferred from the lack of color variation among 
different HII regions (see Section 14.11 and Figure |9j, our SED 
analysis indicates that these star-forming regions form an homo- 
geneous population and that they are significantly different from 



VCC1249 and hence from the two knots belonging to the galaxy 
(C12 and C17). As done previously for VCC1249, by integrating 
over the star formation history, we estimate the stellar masses of 
the studied regions to be between 0.66 x 10 4 and 3.72 X 10 5 M G . 
These values are typically found in star cluster complexes and 
HII regions (Kennicutt et al. 1989). 

6. Discussion 

Our multiwavelength analysis of VCC1249, based on new 
GALEX NUV data, deep optical imaging (including Ha) and 
Keck spectroscopy, reveals that star formation was recently 
quenched in this dwarf and that the only star-forming regions are 
those outside the main body of the galaxy, made of pre-enriched 
gas in the HI cloud. 

These findings are consistent with the conclusions of Sancisi 
et al. (1987), Patterson & Thuan (1992), and Lee et al. (2000) 
that this dwarf has suffered from ram-pressure stripping in the 
hot atmosphere of M49, leading to the ablation of most of its 
original HI gas. Following previous calculations (Sancisi et al. 
1987; Patterson & Thuan 1992) and using our newly estimated 
parameters, we find that VCC1249 is in fact unable to retain its 
gas at any radius because of ram pressure stripping. The radius at 
which ram pressure (Gunn & Gott 1972) becomes efficient can 
be estimated as (Domainko et al. 2006) 

' GM slar M„ as " 

R strip = 0.5R ki - '"" ga J 4 , (1) 
while the stripped mass is 

M strip =M gas {^ + iy X p{-^. (2) 

Here, pHaio = 10~ 3 crrT 3 is the gas density in the halo of M49 
(Fabian 1985) and v = 611 + 32 km s _1 is the relative velocity 
between VCC1249 and M49. In this calculation we furthermore 
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Fig. 14. Summary plot of the photometric points for the HII regions (red), observed optical spectrum if available (green) and their 
best-fit SEDs. Data for CI, C2, C4, C6, and C8 have been modeled using a single burst, while C12 and C17 using an exponentially 
declining SFR. The probability distribution functions for the age parameter of each individual region are shown in the right panels. 



assumed an exponential profile for the stellar and gas compo- 
nents, with 7?o = 14.5 arcsec, the inner scale length computed in 
Section[3]using the g-band luminosity. Adopting M star = 1.20 x 
10 9 M Q , M gas = 8.41 x 10 7 M (since M, 



gas 



■M strip = 1.51 XlO 7 



Mo, Lee et al. 2003), we find that VCC1249 cannot retain its gas 
at any radius, accordingly M strip = M gas . Ram-pressure strip- 
ping can therefore fully deplete the gas reservoir of this dwarf, 
leading to a sudden truncation of its star formatiorQ 

However, as pointed out by Patterson & Thuan (1992) and 
McNamara et al. (1994), ram-pressure stripping is not the only 
mechanism acting on this system because, for instance, it can- 
not explain the presence of the tail and counter-tail of VCC1249 
(Fig. 2b,c). These features instead naturally arise from a tidal in- 
teraction with M49. Patterson & Thuan (1992) estimated that the 
tidal radius, i.e. the distance at which stars in VCC1249 become 
unbound or stripped, is > 20 arcsec, in agreement with the ex- 
cess of light superimposed to the inner exponential profile that 



17 We emphasize that the condition for the ram-pressure stripping is 
satisfied by using the relative velocity between the two galaxies along 
the line of sight. Possible motions in the plane of the sky would make 
ram pressure even more efficient. 



is visible in our deep optical profiles. We conclude therefore that 
both ram pressure and tidal interaction occurred during the in- 
teraction with M49: gravitational tides triggered the tail and the 
counter-tail of VCC1249 and aid ram pressure to remove the HI 
gas by diminishing the overall potential of the dwarf. 

Furthermore, using our SED analysis, we can estimate the 
epoch of the encounter between VCC1249 and M49. Our best- 
fit model for the SED requires a sudden truncation of the star 
formation activity 200 Myr ago because of gas ablation. This 
time is also consistent with an independent estimate for the ab- 
lation time derived by dividing the projected distance between 
the HI cloud and VCC1249 (about 10 kpc) by their relative ra- 
dial velocities (AV = 79 + 30 km s ). Neglecting projection 
effects, this very crude calculation shows that ~ 124 + 47 Myr 
are needed to displace the gas from VCC1249 to the location 
where it is detected in the 21 cm emission. 

Focusing on the outlying HII regions, both the SED fitting 
and the spectroscopy analysis reveal that their stellar popula- 
tions are young and coeval, being born within the last few tens 
of Myr. Furthermore, because all the outlying HII regions have 
an estimated age of less than 30 Myr, which is at least a factor 
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of ten lower than the ablation time, we conclude that these star- 
forming regions were born in situ after the removal of gas. This 
point is strengthened by the observed agreement between the 
oxygen abundance derived spectroscopically for the HII regions 
(12 + log(0/#) = 8.38 for C2 and 12 + log(0/i7) = 8.29forC6) 
with the one for VCC1249 obtained using the mass-metallicity 
relation (12 + log(0/#) = 8.35). 

The presence of a metal-enriched gas tail that extends to 
more than 3 arcmin far from VCC1249 in the direction of M49 
is remarkable also in the context of the metal enrichment of low- 
density gas. Indeed, within dense environments such as rich clus- 
ters, ram-pressure stripping may be an effective way to trans- 
port dust and metals from the ISM of galaxies to the outer intr- 
acluster medium (Boselli & Gavazzi 2006; Cortese et al. 2010). 
Moreover, the encounter of VCC 1 249 with the halo of M49 may 
be the prototype of the interaction between satellite galaxies and 
their centrals and thus this mechanism may be responsible for at 
least part of the metal enrichment seen in the halos of galaxies 
at all redshifts (Tumlinson et al. 201 1). This example shows that 
galactic winds may not be the only processes needed to transport 
metals to large galactocentric distances. 

6. 1 . Witnessing a peculiar interaction 

Witnessing this type of interaction in a cluster must be consid- 
ered fortuitus, a fortiori in the nearby Virgo cluster. Considering 
the morphological types, systems like the VCC1249/M49 pair 
(a giant elliptical galaxy interacting with a dwarf irregular) are 
rare in nearby clusters. In fact, giant elliptical (and cD) galax- 
ies inhabit the center of clusters, while dwarf irregular galaxies 
tend to populate the cluster's outskirts. The majority of dwarf 
galaxies suffer from one or more ram-pressure stripping events 
as they enter the cluster (Bekki 2009), leading to a complete 
removal of the atomic gas and to a consequent transformation 
into dwarf-ellipticals on very short time scales (Boselli et al. 
2008a). Thus, because of a combined effect of a lower veloc- 
ity dispersion and gas density within the cluster, systems like 
VCC1249/M49 would be much more frequent at higher redshift, 
where the population of gas-rich dwarf galaxies has not yet been 
quenched by ram pressure. 

However, looking at low redshift, some similarities with the 
studied system can be found in the famous interacting pair be- 
tween the giant elliptical galaxy (M86) and a giant late-type 
galaxy (NGC 4438) in the Virgo cluster. The highly disturbed 
morphology of NGC 4438 is traditionally (e.g. Boselli et al. 
2005) interpreted as caused by tidal interaction with its compan- 
ion galaxy NGC 4435. Kenney et al. (2008), who discovered an 
extended complex of Ha filaments connecting NGC 4438 with 
M86, instead invoked ram-pressure stripping and tidal (collision) 
interaction with M86. The ram pressure, resulting from the pas- 
sage through the halo of M86, causes the removal of most of the 
ISM from the stellar disk of NGC 4438, while tidal interaction 
produces most morphology disturbances in NGC 4438. As for 
VCC1249, NGC 4438 is very Hi-deficient and what is believed 
to be the remnant of its atomic gas is observed near M86 (Li 
& van Gorkom 2001). Both M86 and M49 have in their prox- 
imity an unusual cloud of HI, left-over of galaxies that passed 
through their halos. However, nothing like the spectacular com- 
plex of Ha filaments connecting NGC 4438 to M86, heated by 
thermal conductivity from the ICM and by turbulent shocks, is 
present between VCC 1249 and M49, where compact Ha emis- 
sion is detected associated with star formation activity taking 
place in compact HII regions. Thus, the presence of an HI cloud 
displaced in the halo of an elliptical galaxy, also combined with 



features visible in the UV-band and Ha, could be the signature 
of a rich-gas galaxy crossing. 

6.2. The complex diffuse structures around M49 

Additional information about the dynamical history of 
VCC 1249 may be gleaned from the tidal debris in the sur- 
rounding field. Early imaging by Patterson & Thuan (1992) and 
McNamara et al (1994) revealed a debris trail ~ 2' in length ex- 
tending northwest from VCC 1249, as well as a short countertail 
to the southwest. More recently, Janowiecki et al (2010) used 
deep imaging to reveal a complex series of accretion structures 
around M49, some possibly associated with VCC1249. With the 
deep NGVS imaging presented here, and the well-determined 
photometric model for M49, we can now search even deeper to 
trace these tidal features further in extent, and search for others. 

To increase signal-to-noise at faint surface brightness 
and highlight any diffuse structures found around M49 and 
VCC 1249, we started with our image of the field after sub- 
tracting the best-fit elliptical isophotal model of M49 (shown in 
Figure[l6ll. We then ran IRAF's OBJMASKS task to identify and 
mask features that are 10cr above a locally defined background 
measured on 1' scales. This effectively masks bright stars, the 
high surface brightness regions of M49's companion galaxies, 
and a myriad of fainter discrete sources on the image. We then 
spatially re-binned the image in 39x39 (7.25"x7.25") pixel bins, 
calculating the median pixel intensity of each bin. The combi- 
nation of masking and medianing effectively rejects the contri- 
butions from small discrete sources and maps the structure of 
diffuse light on ~ 0.5 kpc scales. 

Figure Q3] shows this binned, median-smoothed map. The 
vertical and horizontal banding in the image are artifacts trac- 
ing the edges of the CCDs on the Megacam array. In addition to 
these artifacts, several diffuse features can be seen in the image. 
First, we trace a halo of diffuse light around VCC 1249 out to 
a major axis radius of ~ 2.4' (11 kpc). Beyond this radius, we 
see a stream of light extending from VCC 1249 to the northwest, 
toward the center of M49, coincident with the extended HI and 
UV light. Around M49 itself, there are a number of diffuse fea- 
tures; of particular interest are the shells 5-10' to the northwest 
of M49, originally identified by Janowiecki et al (2010; their 
"inner shell" feature) using deep Schmidt imaging. Because the 
higher spatial resolution of Megacam allows better rejection of 
background contaminants, these shells are more clearly revealed 
in Figure [15] where they are resolved into two features. The in- 
nermost feature extending NW from the center of M49 shows 
a folded fan-like shape, very reminiscent of features that arise 
from radial accretion events (e.g., Hernquist & Quinn 1988). A 
second shell can be seen at slightly larger radius to the west- 
northwest of M49. Other features can be seen in the image as 
well, including more extended shells farther to the northwest and 
to the southeast, and a fan of material extending west-southwest 
of the center of M49. Although the connection between these 
very extended NW and SE shells and VCC 1249 is not clear (as 
Janowiecki et al. 2010 have already pointed out), the better res- 
olution of our NGVS data shows much more detail in the "in- 
ner shell", and clearly links that feature to VCC1249. All these 
features are visible as well in the Schmidt imaging detailed in 
Janowiecki et al (2010); the fact that they are visible in both 
imaging surveys demonstrates that they are not simply artifacts 
generated by instrumental effects such as reflections or scattered 
light. 

We focus here only on the tidal features related to VCC 1249 
itself. In this context, the folded fan-shaped plume to the north- 
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Fig. 15. Left: binned, median-smoothed grayscale map of M49 showing residuals from the ELLIPSE model in the NGVS g-band 
imaging. Although there are artifacts due to the edges of the CCDs and to the bright stars (circular reflection with two sizes: 1 
and 3 arcmin in radius), an extensive series of shells and filaments is apparent in agreement with Janowiecki et al. (2010). Right: 
toy-picture highlighting the diffuse features found in the left-image: the artifacts (green), shells and fan of material described in the 
text (gray) and VCC1249 with its tail pointing toward M49 (blue). 



west of M49 is particularly intriguing. This type of structure 
forms from the radial accretion of the disky companion, as de- 
tailed in Hernquist & Quinn (1988). Indeed, a comparison be- 
tween the structure seen in Figure [l5]and the "t=225" snapshot 
of the HQ88 disk accretion model (their Figure 2) is striking. In 
this scenario, VCC1249 has made (at least) two close passages 
past M49 with the radius of the inner NW shell showing the 
apocentric turning point from the previous passage. This shell 
occurs at a projected radius of 6.4' or 30 kpc from M49's center. 
At this radius, the M49 mass model of Cote et al. (2003) has an 
enclosed mass of 2 x 10 12 Mo, yielding a dynamical time of ~ 
100 Myr. Since VCC1249 is currently found projected 5.6' (26 
kpc) southeast of M49, this is also a rough lower limit on the 
time since the last closest passage to M49. Of course, this limit 
assumes that VCC1249's orbit lies along the plane of the sky and 
the projected separations are true separations, which is likely not 
true. The fact that M49 and VCC1249 have a velocity difference 
of 611 km s~l argues that much of the orbital motion is along 
the line of sight, which significantly increases the inferred time 
since last passage. Given these uncertainties, a scenario where 
VCC1249 is on a fairly radial, bound orbit with last pericenter 
passage occurring a few hundred Myr ago is quite reasonable; 
this orbital timescale is quite similar to the star formation trun- 
cation time we derived earlier (see Section I5.lt and argues for 
a common origirQ In this scenario, tidal stripping of material 
from VCC1249 is ongoing, as is ram-pressure stripping. Indeed, 
because according to this scenario the dwarf galaxy is on its sec- 
ond passage through the dense halo of M49, it likely lost a sig- 
nificant amount of gas on the first passage as well. This first pas- 



18 The timescale since the first passage is much longer than the de- 
rived truncation time. However, the SED analysis can only detect the 
last truncation event, the second passage. 



sage could have stripped the outer, low-density gas in VCC1249, 
leaving the denser gas in the inner regions to be stripped on this 
subsequent passage we observe now. 

6.3. Formation of compact objects 

In addition to the interest for studies of galaxy interactions 
in rich clusters, isolated HII regions born from ram-pressure 
stripped gas provide some insight into the formation of evolved 
compact systems such as globular clusters (GC) and ultra com- 
pact dwarf galaxies (UCD). Found preferentially in the neigh- 
bors of massive elliptical galaxies (Hilker 2011), the origin of 
the latter compact systems is still controversial. Models indicate 
that they can form either by tidal stripping of the diffuse, low sur- 
face brightness disks of nucleated dwarf ellipticals (Bekki et al. 
2003), or might derive from the merging of many young, massive 
(10 6 M ) star clusters, born during a past merger event (Kroupa 
1998; Fellhauer & Kroupa 2002; Kissler-Patig et al. 2006). The 
HII regions formed during the interaction of the (originally) star- 
forming VCC1249 with M49 might be the progenitors of these 
compact systems. They have very compact morphologies and are 
composed of coeval stars formed during a single episode of star 
formation. The stripped gas that gave birth to these HII regions 
is now (dynamically) dissociated from them and thus it cannot 
sustain any more star formation. Given their mass (~ 10 4 M ) 
and their young age (^ 5-10 Myr), it is still unclear whether they 
will survive the "infant mortality" caused by the kinetic energy 
injected by supernova explosions (Lada & Lada 2003). Recent 
models seem to indicate that at least a fraction of them should 
resist and end up in compact evolved systems (Gieles & Bastien 
2008). Containing only 10 4 M , the HII regions found around 
M49 certainly do not have the necessary mass to form UCDs 
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Fig. 16. Grayscale image of M49 showing residuals from the ELLIPSE model that best fits the azimuthally averaged isophotes in the 
NGVS g-band imaging (see Ferrarese et al. 201 1). An extensive series of shells and filaments is apparent. A complex structural was 
also found by Janowiecki et al. (2010); the dashed red lines indicate the regions where these authors found shells and plumes in their 
residual image (here shown for comparison with those in Figure[T5ll. VCC1249 is labeled in cyan, as are VCC1 199 and VCC1 192, 
two compact elliptical galaxies that have likely undergone tidal stripping (e.g., Cote et al. 2010). VCC1205 shows evidence for star 
formation detached from the main body of the galaxy, in the direction of M49. Yellow circles show the position of candidate UCDs 
(having g < 21 and effective radii in the range 10 < R e < 100 pc) identified from the NGVS imaging. At least some of these objects 
may have formed through tidal "threshing" of nucleated dwarf galaxies (e.g., Bekki et al. 2001). 



or massive GCs, but only intermediate-mass compact objects. 
Considering however that at high redshift late-type galaxies had 
gas masses significantly higher than those observed in the local 
universe (e.g. Boselli et al. 2001), the HII regions formed dur- 
ing the interaction of VCC1249 with M49 could be the scaled- 
down version at z — of more massive star clusters formed at 



higher redshift that later transformed into today massive GCs 
and UCDs. Similar star-forming compact structures, but of sig- 
nificantly higher mass, have been indeed observed in the tails of 
some massive spirals in the clusters A 1689 and A2667 at z — 0.2 
by Cortese et al. (2007). 
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6.4. Conclusion 

In conclusion, the interaction between VCC1249 and M49 un- 
derscores the fundamental role played by environment in shap- 
ing the properties of galaxies in cluster cores. In this particular 
case, VCC1249 is undergoing both ram pressure and tidal in- 
teraction: the joint action of the two mechanisms leads to the re- 
moval of the HI gas, while the morphology disturbances are trig- 
gered by the gravitational tides. Our analysis suggests that the 
star formation in VCC1249 was truncated 200 Myr ago, which 
is consistent with the gas ablation time. The HII regions were 
born in situ, within the turbulent, pre-enriched gas that was re- 
moved by the interaction. 

6.5. A harsh environment: the Virgo B subcluster 

We conclude with a panoramic view of the immediate neigh- 
borhood of VCC1249 — an environment that is, of course, 
dominated by M49, the brightest member of the Virgo cluster 
and the central galaxy in Virgo's B subcluster (e.g., Binggeli, 
Popescu & Tammann 1993). Our analysis of VCC1249 adds 
to the growing body of evidence that interactions, mergers and 
stripping have had a profound affect on the galaxy population 
in this high-density environment. Figure [16] illustrates some of 
this evidence. In addition to VCC1249, this region also contains 
two examples of the rare class of "compact elliptical" galaxies, 
whose origin is almost certainly related to strong tidal stripping 
of initially more massive galaxies (e.g., Faber 1973; Bekki et al. 
2001b; Chilingarian et al. 2009; Cote 2010; Huxor et al. 201 1). 
Furthermore, as shown in the residual image in Figures [TBI and 
[TBI a complex series of shells, plumes and streams surrounds 
M49 — prima facie evidence for past accretions and interactions 
(Janiowiecki et al. 2010; Ferrarese et al. 2011). Fitting of PSF- 
convolved models to sources in the NGVS images also reveals a 
large number of UCDs candidates, which are shown as the yel- 
low circles in Figure[16](see also Hasegan et al. 2005). Although 
there may well be multiple formation channels for UCDs, at least 
some of these objects could have formed through tidal "thresh- 
ing" of nucleated dwarf galaxies, a leading UCD formation 
mechanism (Bekki et al. 2001a). Finally, approximately 11' to 
the south of M49 lies VCC1205, which is classified as ScIII-pec 
by Binggeli et al. (1985). Despite its high relative velocity with 
respect to M49 (Av = VM49 ~ Vvcct205 = 1001-2341 = 1340 km 
s i), the galaxy shows a minor HI deficiency, while its optical 
and NUV morphology shows evidence for star-forming regions 
detached from the main body of the galaxy in the direction of 
M49 (Figure [T7li. This is perhaps another indication of the harsh 
environment existing in the core of the Virgo subcluster B. In 
any case, there is mounting evidence that the central regions of 
rich clusters — and of the Virgo B subcluster in particular — are 
highly dynamic environments in which interactions influence the 
structure of both the central galaxies and their surrounding satel- 
lites. 
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